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Abstract 

Using the n-particle periodic Toda lattice and the relativistic generalization due to 
Ruijsenaars of the elliptic Calogero-Moser system as examples, we revise the basic 
properties of the Backlund transformations (BT's) from the Hamiltonian point of view. 
The analogy between BT and Baxter's quantum Q-operator pointed out by Pasquier 
and Gaudin is exploited to produce a conjugated variable fi for the parameter A of 
the BT B\ such that fi belongs to the spectrum of the Lax operator L(X). As a 
consequence, the generating function of the composition B\ 1 o . . . o B\ n of n BT's gives 
rise also to another canonical transformation separating variables for the model. For 
the Toda lattice the dual BT parametrized by jd is introduced. 
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1. Introduction 



Backlund transformations (BT's) are an important tool in the theory of integrable sys- 
tems |J . Most frequently, they are understood as special mappings between solutions 
of nonlinear evolution equations. The Hamiltonian properties of BT's, as canonical 
transformations, are studied less well. The recent developments in the quantum inte- 
grable theories 0, |3| , discrete-time dynamics ||, ||, |J and separation of variables |7|, |[ 
suggest, however, that the Hamiltonian aspect of BT's deserves more attention. 

The aim of the present paper is to revise the concept of BT's from the Hamiltonian 
point of view and to point out some new properties of BT's. We restrict our attention 
to the finite-dimensional integrable systems and illustrate our general remarks on the 
example of the periodic Toda lattice and the elliptic Ruijsenaars model. When elab- 
orating our approach to BT's, we have benefited greatly from the works of Pasquier 
and Gaudin 0, where a fundamental relationship between BT and Baxter's quantum 
Q-operator was discovered, and of Veselov Q, who gave us the adequate mathematical 
language to speak about integrable mappings. 

In the section 2 the main properties of Backlund transformations for Liouville in- 
tegrable systems are enlisted and a new property of spectrality is introduced. The 
meaning of spectrality is elucidated by making the comparison with the Baxter's quan- 
tum Q-operator. It is shown that spectrality of BT provides an effective solution to the 
problem of separation of variables. In two subsequent sections we illustrate the new 
property of BT's for two families of integrable many-body systems. The concluding 
section 5 contains a summary and a discussion. 



2. Spectrality and separation of variables 

Suppose an integrable system with n degrees of freedom is described in terms of the 
canonical Darboux variables X = {X{\™ =1 and x = {xj}" =1 , with the Poisson brackets: 

{Xi,Xj} = {xi,Xj} = 0, {X h Xj} = Sij, (2.1) 

and functionally independent commuting Hamiltonians Hi = Hi(X,x) 

{H l ,H j } = 0, i,j = l,...,n. (2.2) 

For our purposes it is convenient to think of a BT as a canonical transformation B\ 
from the canonical variables (X, x) to the canonical variables (Y, y) . It is important 
that B\ depends on a complex parameter A. We shall suppose that B x can be described 
via the generating function F\(y\ x) such that 

X- 9Fx Y- 9Fx (2 3) 

Xi ~ ^? Yl ~ ~^7- (2 ' 3) 

The list of properties defining a BT usually includes: 
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• Canonicity. See above. 

• Invariance of Hamiltonians. 

H i (X,x) = H i (Y,y), i = l,...,n. (2.4) 

• Commutativity. 

B Xl oB X2 = B X2 oB Xl (2.5) 
where o means composition of canonical transformations. 

In case of the algebraically integrable systems |J one more property can be added 
to the list: 

• Algebraicity. The equations ( |2.3|) describing B x are supposed to be algebraic with 
respect to X, Y and properly chosen functions of x and y (say, exponential or 
elliptic). 

In the present paper, however, we concentrate on the analytic properties of BT's 
and ignore their algebraic and algebro-geometric aspects. 

It is important to make clear distinction between the notion of BT and the close 
notions of integrable canonical mapping or integrable discrete-time dynamics. The 
latter two are defined by the properties of canonicity and invariance only, the parameter 
A being disregarded. The term 'discrete-time dynamics' refers usually to the case 
when the canonical transformation degenerates, in a certain limit, into an infinitesimal 
generator {H, •} of a continuous Hamiltonian flow. Existence of the parameter A is 
crucial for our definition of BT and enriches it with new properties. 

Though the commutativity of BT's is traditionally proved as an independent prop- 
erty, in fact it follows from the canonicity and the invariance of Hamiltonians. Indeed, 
as shown in ||, any integrable canonical mapping acts on the Liouville torus as a 
shift (or a collection of shifts, in case of multivalued mappings) of the angle variables 
<fi — > tpi + bi(X). The commutativity is then obvious. 

The theory of BT's acquires a new aspect if the integrable system in question is 
solvable via Inverse Scattering (or Inverse Spectral Transform) method. Suppose that 
the commuting Hamiltonians Hi can be obtained as the coefficients of the characteristic 
polynomial 

W(u, v; {H { }) = det(v - L{u)) (2.6) 

of a matrix L(u) = L(u; X, x) (Lax operator) depending on X, x and a complex pa- 
rameter u. Note that the invariance of Hi under B x is equivalent then to the invariance 
of the spectrum of L(u), that is there exists an invertible matrix M(u) such that 

M(u)L(u;X,x)=L(u;Y,y)M(u), Vu G C. (2.7) 

The properties of BT's enlisted above are well known. Now we are going to add to 
the list a new property which is the main contribution of the present paper. 
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• Spectrality. Let \i be defined as the variable conjugated to A: 

We shall say that the BT B\ is associated to the Lax operator L(u) if for some 
function /(/i) the pair (A, /(/x)) lies on the spectral curve of the Lax matrix 

W(A, /(//); = det(/(/i) - L(A)) = 0. (2.9) 

This spectrality property of BT seems to be new, at least we failed to find it in 
the literature. We have verified it for the Toda lattice and the elliptic Ruijsenaars 
model for which /(/i) = (see sections 3 and 4). It seems plausible, however, that 
spectrality is the property shared by BT's for a much larger class of models. 

The meaning of the equality ( |2.9| ) becomes clear if we turn to the quantum case. 
In the pioneering paper by Pasquier and Gaudin 0, based on the earlier treatment of 



the classical Toda lattice by Gaudin [10J, a remarkable connection has been established 
between the classical BT B\ for the Toda lattice and the famous Baxter's Q-operator 
|fnj . Pasquier and Gaudin have constructed certain integral operator Q\ 

Qx : *(x) -> J dx Q x (y; x)^(x) (2.10) 

(here and below dx = dxi A ... A dx n etc) whose properties parallel those of the 
classical BT B\. In the quantum case the canonical transformation is replaced with 
the similarity transformation 

Y i = Q x X l Ql\ y i = Qxx l Q- x 1 , (2.11) 

where the hat A distinguishes the quantum operators from their classical counterparts. 
The correspondence between the kernel Q\(y; x) of Q\ and the generating function 
F\(y; x) of B\ is given by the semiclassical relation 



Qx(y) x) ~ exp (~F x {y; ^)) , h ^ 0. (2.12) 



After publication of [[| the Q-operators have been found for a number of other 
quantum integrable models 0. 

The properties of Q\ such as the invariance of the Hamiltonians 

[Q Xi Hi] = (2.13) 

and the commutativity 

[Qm,Qx 2 ]=0 (2.14) 

reproduce the respective properties (|2.4j ) and ( |2 . 5| ) of B\. The most interesting property 
of Q\, however, is that its eigenvalues 0(A) on the joint eigenvectors ^ u of Hi and Q\ 
labelled with the quantum numbers v 

Qx* v = 4> v {\)* v (2.15) 
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satisfy the separation equation, which is a certain differential or difference equation 



w(\,-ih^,{h i }U v (\)=0 (2.16) 



containing the eigenvalues hi of Hi. In the classical limit the equation ( 2.1 6|) goes over 
into the spectrality equation ( [2.9|) . 



An important application of the spectrality property of BT is that to the problem 
of separation of variables |7], H]- Again, it is instructive to start with the quantum 
case. A separating operator K is, by definition, an operator, transforming the joint 
eigenfunctions ^ v of Hi into the product 

Ky„ = c„f[cj) u (^) (2-17) 

i=l 

of separated functions 4> U {X) of one variable A satisfying the separation equation ( ^.16| ). 
Since the coefficients c v in ( |2.17|) can be chosen arbitrarily, abstractly speaking, there 



exist infinitely many separating operators K. The difficult problem, however, is to find 
the ones which can be described as integral operators with explicitely given kernels. 

Knowing a Q-operator gives one an immediate opportunity to construct plenty of 
separating operators. Indeed, consider the operator product Q\ x ,.,\ n = Qx 1 ---Qx n 
having the kernel Q\ 1 ...\ n (y; x) and for any function p{y) introduce the operator 

K p : *(x) - JdxJ dyp(y)Q Xl ... Xn (y;x)^(x). (2.18) 

It is obvious from ( |2.15|) that K p is a separating operator, the coefficients c v being 

c v = J dyp(y)V u (y). (2.19) 

Since the eigenfunctions ^ v {y) form a basis in the corresponding Hilbert space, the 
formula ( |2.19| ) provides a one-to-one correspondence between reasonably chosen classes 
of c u and p{y). Therefore, arguably, the formula ( [2.18| ) describes all possible separating 
operators. Their kernels K p (X; x) are given explicitely as multiple integrals 

K p (X;x) = JdyJd^...Jd^ 

xp(y)QxAy^ (1) )Qx 2 (e ) ;t {2) )---Qx^ {n - 1) ^)- (2.20) 

It is a straightforward task to present the classical analog of the above argument. 
Consider the composition Bx 1 ...\ n = B\ 1 o . . . o B\ n of Backlund transformations and 
the corresponding generating function i 7 A 1 ...A n (?/; x). Let us switch now the roles of ?/s 
and A's treating A's as dynamical variables and y's as parameters. Then F\ 1 ^\ n (y; x) 
becomes the generating function of the n-parametric canonical transformation K y from 
(X, x) to (/i, A) given by the equations 

Xi - ~d^r> * = — d\r • (2 - 21) 
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It follows directly from ( |2.9| ) that the pairs (Aj,/ij) satisfy the separation equations 



W(X i ,f(/jL i y,iH j }) = Q (2.22) 

which constitutes exactly the definition of the separating canonical transformation in 
the classical case 0. 

The above construction corresponds in the quantum case to setting p(y) = 8(yi — 
yi) . . . 5{y n — y n ) where y~i are some constants. It remains an open question what could 
be the classical analog of the formula ( [2.18|) for generic p{y)- 



As the last general remark before passing to the examples, we would like to stress 
that for the finite-dimensional systems the composition of n BT's with n being the 
number of degrees of freedom is a sort of 'universal' BT in the sense that any other 
canonical transformation preserving the Hamiltonians Hi must be expressible in terms 
of Bx L ...\ n - To observe it one can use again the fact that in the angle coordinates B\ 
acts as a shift ipi — > tpi + bi(\). For generic bi(\) the sum bi(Xi) + . . . + 6j(A n ) must then 
cover the n-dimensional Liouville torus which results in the universality of £?Ai...A„- 



3. Periodic Toda lattice 



Our first example is the periodic Toda lattice [|12j , |13|| for which there exist two alterna- 
tive Lax operators associated, as we shall show, with two different BT's. The standard 
and quite well studied BT [|TJ, [T2|, [Tl [ITJ which we denote here B\ is associated, in the 



sense defined in the previous section, to the 2x2 Lax matrix (or, monodromy matrix 
[151) L(u; X, x) defined as the product of local L-operators 



L{u)=£ n {u)...£ 2 {u)£ 1 {u), (3.1) 

ti(u) = tifaXi,*) = ( U + x Xl ~f ) . (3.2) 

The characteristic polynomial of L(u) is quadratic in v 

W{u, v) = det(> - L(u)) = v 2 - t(u)v + 1, (3.3) 

and the commuting Hamiltonians Hi are obtained from the expansion of the only non- 
trivial spectral invariant t(u) = tvL(u) 

t(u) = u n + H x u n - X + ... + H n . (3.4) 

In particular, 

1 n /I \ 

-Hl-H 2 = Y { {- 2 X 2 + ^-^ (3.5) 

is the standard periodic Toda Hamiltonian (in this section we use the periodicity con- 
vention % + n = i for the indices i). 
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The Backhand transformation B\ is obtained from the generating function 

n 

Fx(y; x) = ]T(e^ - e^ 1 "*' - X(x l - Vi )) (3.6) 



and, according to ( |2.3| ), is implicitely described by the equations 

Xi = e Xi ~ Vi + e yi+1 - Xi -A, Yi = e Xi ~ Vi + e" 4-3 "- 1 - A. (3.7) 

The characteristic properties of the BT are verified easily. The invariance of the 
Hamiltonians can be established using the equality [ |ITJ[ 

M i+ i(u, \)li(u; X u — £i(u; Y h y^M^u, A), (3.8) 

where 

/ 1 —QVi \ 

Mi(u,X) = Mi(u,\;xi-x,yi) = I e - Xi _, X _ u _ e y i - Xi - 1 j ( 3 - 9 ) 

which one can verify directly using the equations ( |3.7| ). Due to the periodic boundary 
conditions, the local gauge transformation ( |3.8|) results in the spectrum-preserving 
similarity transformation 

Mi(-u, X)L{u; X, x) = L{u;Y,y)M 1 (u,X) (3.10) 

for L(u) which proves the invariance (|2.4|) of the Hamiltonians. 



The direct proof of the commutativity (2J3) of the BT's can be found in |T3|, [12], |T|. 

To prove the spectrality equality (|2.9|) which in this case takes the form det(e _/i — 
L(A)) = we shall apply a modified version of the argument used in [H], 0] for the 
quantum case. Note, first, that in our case 

dF x n 

v = --Q^=J2( x i-yi)' ( 3 - n ) 

as follows from ( |3.6j ) and ( |2.8| ). It suffices then to show that e~ M is an eigenvalue of the 
matrix L(\). We shall construct explicitely the corresponding eigenvector u\\ 

L(X;X 1 x)u 1 = e-' i u 1 . (3.12) 

From ( |3~^| ) it follows that det(Mj(u, A)) = A — u. It is easy to see that for u = A 
the matrix Mj(A, A) has the unique, up to a scalar factor, null- vector 

^=^ 6 ^, Mi(\,\)ui = 0. (3.13) 

Using the identity ( |3.1(J| ) we conclude that 

Mi(A, A)L(A; X, x)u x = (3.14) 

which, combined with the uniqueness of the null- vector a>i of Mi, implies that lo\ is 
an eigenvector of L(A; X, x). To determine the corresponding eigenvalue, we apply the 
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same argument to the identity ( |3.8| ) obtaining the equality M i+ i(X; X)£i(X; X{, Xi)u>i = 
from which it follows that £i(X; X{, Xi)uJi ~ u>i + i. The direct calculation shows that 

£i(X; Xi, Xi)uJi = e yi ~ Xi uj i+1 . (3.15) 

It remains only to use the formulae (|3.1| ) and ( |3.11| ) to arrive finally at ( p,12|) . Actually, 
we could skip the discussion of null-vectors of M; and to derive ( |3.12|) directly from 
( |3.15| ). In more complicated situations, however, it may be easier to find u as the 
null- vector of M and then to determine the corresponding eigenvalue of L(X). 

Note that the vectors u>i are the classical counterparts of Baxter's [[TTJ] vacuum 
vectors. 

Let us examine now the alternative Lax operator [12|, [TjJ given by the n x n matrix 
C{y; X, x) with the components 

C jk (v; X, x) = -Xfa + v- 1/n e x >- Xk 5 hk+1 + v l/n 5 ]+hk . (3.16) 

The duality between the Lax operators C(v) and L(u) is expressed in the switching 
the roles of the parameters u and v. The characteristic polynomial W(v, u) = det(-u — 
C{v)) of the Lax operator (|3.16| ) produces the same Hamiltonians Hi and the same 
spectral curve as W(u,v), as follows from the identity 

det(v - L(u)) = -vdet(u - C(v)). (3.17) 

For other examples of the similar duality, see JM. 

The swapping of u and v corresponds to switching the roles of the parameters A and 
/j, in the BT. For the new Backlund transformation £> M associated with the Lax operator 
C{v) the formulae ( |3.6| ), ( p.7|) and ( |3.11|) remain the same but their interpretation 
changes. The BT is parametrized now by the parameter \i which becomes a numerical 
constant. The equality ( 3.11|) is reinterpreted now as a constraint on the variables Xi 
and yi. The parameter A is reinterpreted, respectively, as the Lagrange multiplier for 
the constraint ( |3.11| ) and becomes a dynamical variable which can be defined from the 
equations ( |3.7| ). 

The characteristic properties of BT are verified for in very much the same 
manner like for B\. The invariance of the Hamiltonians follows from the invariance of 
the spectrum of C(v) which, in turn, follows from the easily verified identity 

M(v)C(v; X, x) = C(v; Y, y)M(v), (3.18) 

with the matrix Ai(v) = Ai(v; x, y) given by its components 

M jk (v) = -8 jk + v- l ' n ^- x *5 hk+l . (3.19) 

The commutativity, as shown in section 2 follows from the canonicity and the in- 
variance. 

To prove the spectrality equality det(A — £(e _/i )) = 0, it suffices, similarly to the 
case of the Lax operator L(u), to present the eigenvector Q of the matrix £(e _/i ) 
corresponding to the eigenvalue A: 

C(e-»)Q = XVI. (3.20) 
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Again, Q can be determined as the null- vector of A4(e fl ): 



Ml 



)Q = 0. 



(3.21) 



Note that the uniqueness of Q follows from the easily verified identity det(z — Ai(v )) = 
(z + l) n — v~ l e~ fl which implies that the spectrum of .M(e -M ) consists of n non- 
degenerate eigenvalues, the being one of them. From (|3.21| ) one easily derives the 
recurrence relation for the components of Q 



= fij-i exp I yj 



+ * 

n 



(3.22) 



which determine f2 up to a constant factor. It remains to verify the identity ( |3.20 ) 
which can be done by a direct calculation using the expressions ( p. 16 ) for the matrix 
C(v), (fm]) for n and Q for X t . 



4. Elliptic Ruijsenaars model 

Our second example is the relativistic generalization due to Ruijsenaars [T7J] of the el- 



liptic Calogero-Moser [18| many-body system. For the non- relativistic Calogero-Moser 



system a BT was found in |I9 1 . In p] a discrete-time dynamics was constructed for the 
elliptic Ruijsenaars model. As we show below, the discrete-time evolution transforma- 
tion found in || has all the properties of a BT if the parameter p in || is specified in 
a proper way. 

We use here the notations of [||] with few exceptions: our parameter £ equals to 
—A from P, and our e _A corresponds to p from As in the case of the Toda lattice, 
there exist two dual BT's: B x and £y The standard Lax operator for the Ruijsenaars 
model, as shown below, is associated with B^. Since the dual Lax operator is so far 
unknown, we describe here only the transformation B^. 

Following R M, we introduce the Lax operator C(v;X,x) for the n-particle (A n _i 
type) Ruijsenaars system as the n x n matrix with the entries 

{v) = _ eX MQa(y + ^ - ^ - n ^ ; -x fc + (41) 

a{v)a{Xi - Xj - Q a[Xi - x k ) 

where a(x) is the Weierstrass sigma function and £ is a constant. 
The commuting Hamiltonians 

JC(l,..,n} \jeJ J J6J U \- L j X k) 

\J\=i V ' fce{l,-..,n}\J 

are generated from the characteristic polynomial of the matrix C(v) (|4.1|) 

det(C(v) -u) = Yi-uV^Hj y - ^ (4.3) 
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where we assume Hq = 1. 

The Backlund transformation is given by the equations 

e* = e"* f[ ^ = ~ + (4-5) 
where A is considered as the Lagrange multiplier corresponding to the constraint 

n 

£* = «£ + XX^fc -Vk)- ( 4 -6) 
fc=i 



Note here that the variable A in formulas (|4.4|) - (|4~o|) , describing the dicrete-time dy- 
namics, appeared already as p in ||, but the conjugated variable \i did not. Notice 
also that A was treated in pj as an extra parameter, not as a Lagrange multiplier 
corresponding to a constraint. 

The generating function of the canonical transformation £> M is expressed in terms 
of the function x 

S{x) = J In a{y)dy (4.7) 

as follows: 

n 

F\(y;x) = -^^2{.x i -y i + ^)+^2( l S(x i -x j -^)-S{xi-x j + ^)) 

i=X i<j 
n 

+ J2(S(x i -y j + 0-S(x i -y j )). (4.8) 

The verification of the characteristic properties of BT for £> M proceeds in the same 
way as in the case of the Toda lattice. 

The invariance of the Hamiltonians Hi follows from the identity (see || for the 
proof) 

M(v)C(v;X,x) — C(v;Y,y)M(v) (4.9) 
where the matrix A4(v) = Ai(v; x,y) is defined as 

a{Vi - xj - f ) a(yi - y k ) y a(x k - 

The commutativity, as usual, is a consequence of canonicity and invariance (see 
section 2). 

To prove the spectrality equality which takes the form det(e _A — C(ji)) = it is 
sufficient, like in the case of the Toda lattice, to find the eigenvector Q of the matrix 
C(fj) corresponding to the eigenvalue e~ A . Let us show that, up to a constant multiplier, 
the components of the eigenvector Q are 



fij = ^ ; — . (4.11) 

n o-(xi - x k ) 
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The equality 



£(/z)fi = e~ x tt 



or 



3=1 



(4.12) 



(4.13) 



after the substitutions (|4.1|) for Cij, ( |4.6|) for /x and (|4.4|) for e Xi is reduced to the 
following identity for sigma functions 

£ <t(j* + x t - x 3 - Oil -Vk + OU ( ~ ^ ~, =a^)Y[a(x t -y k ). (4. 14) 



3=1 



k=l 



a[Xj - x k 



k=l 



Due to the symmetry, it is sufficient to prove ( |4.14|) only for i = 1. Let i = 1 and 
n > 2. Consider the both sides of the equality Q4.14Q as functions of x n . It is easy 
to see that they are holomorphic in x n (the apparent poles in the left-hand-side being 
cancelled) and have the same quasiperiodicity properties. From the basic properties 
of sigma functions |20[ it follows that it is sufficient to verify the equality of LHS and 
RHS only in one arbitrary point x n = x with the only condition x ^ fi — x n . Choosing 
x = y n — £ we observe that ( J4.14 ) is reduced to the similar identity of order n — 1 . The 
proof follows then by induction in n since the case n — 1 is trivial. 

As in the section 3, the vector Q ( [4.1 1| ) is again the null- vector of the matrix Ai(fi), 

i.e. 

M(n)Q = 0. (4.15) 
The corresponding identity for sigma functions 



2^ a W + Vi~ Xj - ; — = 



3=1 



Uk^j o{xj - x k ) 



(4.16) 



follows from the identity 



E |! L f 3 Zk \ =0 if ±{zk-x k ) 

j=1 Uk^j <?(xj ~ x k) 



(4.17) 



k=l 



(Cf. 



p. 451) when one substitutes z k = y k — £ for k ^ i and % = fi + — £. 



5. Discussion 



We have studied three new aspects of Backlund transformations. Those are spectral- 
ity dual BT's and application of BT's to the problem of separation of variables. As 
demonstrated in the section 2, the composition of n BT's, being an 'universal' (n- 
parametric) BT, provides a separation of variables which has n arbitrary parameters, 
and thereby defines an 'universal' (n-parametric) family of separating transformations. 
The connection between the 'universal' BT and the 'universal' SoV is intriguing and 
has yet to be studied in detail. 
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Though we have discussed in the present paper only the classical case, our pri- 
mary motivation comes from the quantum case. The main problem in the quantum 
case is to construct Baxter's Q-operator which is a quantum analog of the Backhand 
transformation. For the trigonometric case of the Ruijsenaars system, i.e. for the case 
of multivariable (y4 n _!-type) Macdonald polynomials, we have succeeded to describe 
explicitely such a quantum analog of the transformation introduced in the section 
4. The results will be reported elsewhere. 
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